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We present a measurement of CP-violating asymmetries in fully reconstructed B°^D^ ± n^ : 
decays in approximately 88 million T(4S) — > BB decays collected with the BABAR detector at 
the PEP-II asymmetric-energy B factory at SLAC. From a time-dependent maximum likelihood 
fit we obtain for the CP-violating parameters: a — —0.022 ± 0.038 (stat.) ± 0.020 (syst.), a* = 
-0.068±0.038 (stat.) ±0.020 (syst.),ci ep = +0.025±0.068 (stat.) ±0.033 (syst.), and cj* ep = ±0.031± 
0.070 (stat.) ± 0.033 (syst.). Using other measurements and theoretical assumptions we interpret 
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the results in terms of the angles of the Cabibbo-Kobayashi-Maskawa unitarity triangle, and find 
| sin(2/3 + 7)| > 0.69 at 68% confidence level. We exclude the hypothesis of no CP violation (sin(2/3 + 
7) = 0) at 83% confidence level. 

PACS numbers: 12.15.Hh, 11.30.Er, 13.25.Hw 



In the Standard Model, CP violation in the weak inter- 
actions between quarks manifests itself as a non-zero area 
of the Cabibbo-Kobayashi-Maskawa (CKM) unitarity tri- 
angle . While it is sufficient to measure one of its angles 
a, ft, or 7 to be different from or 180° to demonstrate 
the existence of CP violation, the unitarity triangle needs 
to be overconstrained with different measurements to test 
the CKM mechanism. Measurements of (3 free from the- 
oretical uncertainties exist 0, Q , but there are no such 
measurements of a and 7. This letter reports the mea- 
surement of CP-violating asymmetries in B°—>D^ ± t: t 
decays in T(4S) — > BB decays and its interpretation 
in terms of constraints on | sin(2/3 + 7)! 0,0. 

The time evolution of B°^D^ ± t: t decays is sensi- 
tive to 7 because of the interference between the CKM- 
favored decay B°— >Z?(*) + 7r~ , whose amplitude is pro- 
portional to the CKM matrix elements V cb V* d , and 
the doubly-CKM-suppressed decay B ^D^ + tt~ , whose 
amplitude is proportional to V cd V* b . The relative weak 
phase between the two amplitudes is 7, which, when com- 
bined with B°B° mixing, yields a weak phase difference 
of 2/3 + 7 between the interfering amplitudes. 

The decay rate distribution for B°^D ± n T decays is 



At) 



-|At|/T 

4r 



(1) 



[1 T S c sin(Am d At) =p r)C cos(Am d At)] , 



where r is the B° lifetime, neglecting the decay width 
difference, Am^ is the B°B° mixing frequency, and 
At = t rec — ttag is the time of the B a — ► D ± tt t decay 
(B lec ) relative to the decay of the other B (Stag). In 
this equation the upper (lower) sign refers to the flavor 
of as B° (B°), while rj = +1 (-1) and ( = + (-) for 
the final state D~ir + (D + tt~). In the Standard Model, 
the <S* and C parameters can be expressed as 



S± 



2Im(A±) 
1 + |A±| 2 



and C 



1 



(2) 



where A± = r 1 e~H 2 P+7T Here S is the relative strong 
phase and r is the magnitude of the ratio of the sup- 
pressed and the favored amplitudes. The same equations 
apply for B°^D* ± it t decays, with r and S replaced by 
the parameters r* and <5*, respectively 0. 

The analysis strategy is similar to that of the time- 
dependent mixing measurement performed at BABAR || . 
To identify the flavor of Stag, each event is assigned by 
a neural network to one of four hierarchical, mutually 
exclusive tagging categories: one lepton and two kaon 
categories based on the charges of identified leptons and 



kaons, and a fourth category for remaining events. The 
effective tagging efficiency is (28.1 ± 0.7)% The time 
difference At is calculated from the separation along the 
beam collision axis, Az, between the B rec and Stag decay 
vertices. We determine the B lec vertex from its charged 
tracks. The S tag decay vertex is obtained by fitting 
tracks that do not belong to B lcc , imposing constraints 
from the B rcc momentum and the beam-spot location. 
The At resolution is approximately 1.1 ps. 

The expected CP asymmetry in these decays is small 
(r« « \V: b VjV: d V cb \ « 0.02), and therefore this mea- 
surement is sensitive to the interference between the 
b — > u and b — > c amplitudes in the decay of S tag . To 
account for this effect we use a parametrization different 
from Eq. [21 which is described in Ref. || and summa- 
rized here. For each tagging category (i) the interference 
is parametrized in terms of the effective parameters r[ 
and 5[. Neglecting terms of order r'*' 2 and r,- 2 , for each 
tagging category the At distribution becomes 



,-\M\/-r 
At 



[1 T (a W T r]bi - r/c. 



sin( Am d At) =p V cos(Am d At)] , (3) 



where, in the Standard Model, 

a W = 2r ( * ) sin(2/3 + 7)cos(S ( * ) , 
bi = 2r-sin(2/3 + 7)cos(5,', 

c ( * ] = 2cos(2/5 + 7)(r ( * ) sin<5 ( * ) -r-sin^). (4) 

Semileptonic B decays do not have a doubly-CKM- 
suppressed amplitude contribution, and hence r[ cp = 0. 
Given that we have two B decay modes and four tag- 
ging categories, we use two a parameters (one for each 
final state), three b parameters (one for each non- lepton 
tagging category), and eight c parameters (one for each 
combination of tagging category and final state). Re- 
sults are quoted only for the four parameters a*-*- 1 and 
c[*p, which are independent of the unknowns r\ and 5[. 
The other parameters are allowed to float in the fit, but, 
as they depend on r[ and 6^, they do not contribute to 
the interpretation of the result in terms of sin(2/3 + 7). 

This measurement is based on 88 million T(AS) — > 
BB decays, corresponding to an integrated luminosity 
of 82 fb^ 1 , collected with the BABAR detector 1Q| at the 
PEP-II asymmetric-energy B factory at SLAC. We use a 
Monte Carlo simulation of the BABAR detector based on 
GEANT4 11] to validate the analysis procedure and to 
estimate some of the backgrounds. 

The event selection and the reconstruction of B° — > 
£)(*)± 7r T candidates are detailed in Ref. 0. Signal 
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FIG. 1: Distributions of J71.es in the A_E signal region for 
events with tagging information in the B° — » D ± 7r =F (left 
plot) and the B° — > D**-^ sample (right plot). 



and background are discriminated by two kinematic 
variables: the beam-energy substituted mass, toes = 

VW2) 2 



p* B , and the difference between the B can- 
didate's measured energy and the beam energy, AE = 
Eg — (y / i/2), where Eg (p* B ) is the energy (momentum) 
of the B candidate in the e + e~ center-of-mass frame, 
and y/s is the total center-of-mass energy. The signal 
region is defined as \AE\ < 3<t, where the resolution a 
is mode-dependent and approximately 20 MeV, as deter- 
mined from data. Figure[l]shows the toes distribution for 
candidates in the AE signal region. The toes distribu- 
tion is fit with the sum of a threshold function 12], which 
accounts for the background from random combinations 
of tracks, and a Gaussian distribution with a fitted width 
of about 2.5MeV/c 2 describing the signal. After tagging, 
the Gaussian yield is 5207 ± 87 and 4746 ± 78 events 
for the B° — > D^n 11 and B° — > D* ±r K T sample respec- 
tively, with corresponding purities of (84.9 ± 0.5)% and 
(94.4±0.4)% in a ±3<r region around the nominal B mass. 
Backgrounds from B a decays that peak in the toes sig- 
nal region were estimated with Monte Carlo simulation 
to constitute (0.21 ± 0.06)% and (0.13 ± 0.05)% of the 
B° —>■ D ± n T and B a — > D* ± tt t yields, respectively. For 
backgrounds from B + decays, the corresponding figures 
are (0.93 ± 0.23)% and (0.93 ± 0.10)%. 

An unbinned maximum-likelihood fit is performed on 
the selected B candidates using the At distribution in 
Eq. [3] convolved with a resolution function composed 
of three Gaussian distributions. Incorrect tagging di- 
lutes the parameters a^*\ cf\ and the coefficient of 
cos(ArridAt) by a factor Di = 1 — 2wi 0, 0, where 
u>i is the mistag fraction. The resolution function and 
the parameters associated with flavor tagging are deter- 
mined from the data and are consistent with previous 
BABAR analyses 01 • The combinatorial background is 
parametrized as the sum of a component with zero life- 
time and one with an effective lifetime fixed to the value 
obtained from simulation. The fraction of each compo- 
nent and the At resolution parameters are left free in 
the fit to the data. The background coming from B^ 
mesons is modeled with an exponential decay with the 
B^ lifetime, and its size is fixed to the value predicted 



by simulation. The background from B° mesons is ne- 
glected in the nominal fit, but is considered in evaluating 
the systematic uncertainties. 

The results from the fit to the data are 

a = -0.022 ± 0.038 (stat.) ± 0.020 (syst.) , 

a* = -0.068 ±0.038 (stat.) ± 0.020 (syst.), 

ci C p = +0.025 ±0.068 (stat.) ± 0.033 (syst.), 

ci* cp = +0.031 ± 0.070 (stat.) ± 0.033 (syst.) . (5) 

All other fitted b and c parameters are consistent with 
zero. Figure[5]shows the fitted At distributions for events 
from the lepton tagging category, which has the lowest 
level of background and mistag probability. 

The systematic uncertainties on the parameters in 
Eq. [S] has been calculated in a manner similar to that 
used in Ref. A small bias in the At measurement 
could result in a bias on the c parameters in Eq. [3J For 
instance, a realistic At bias of 0.024 ps results in a shift in 
c fe P °f 0.002. We are immune from this effect because we 
fit for tagging category dependent biases in the resolution 
function directly on data. Nonetheless, the impact of a 
possible mismeasurement of At has been estimated by 
varying the assumptions on the resolution function, the 
position of the beam-spot, the absolute z scale, the inter- 
nal alignment of the vertex detector, and quality criteria 
on the reconstructed vertex. The corresponding error on 



0.015, while that on cW is a c = 0.026. The 



ttW IS O a 

systematic uncertainties on the fit technique [a a = 0.013, 
cr c = 0.020) include the upper limit on the fit bias esti- 
mated from samples of fully simulated events, the uncer- 
tainty on the B° lifetime and Ato^ |l3| . and the impact 
of neglecting higher order terms in A*^ or r[ in Eq. [3] 
As a cross-check, we performed the same fits on sam- 
ples of 18233 B~ -> DO) ?!-" and 1740 B° -> J/ip K*° 
candidates, where we find no significant CP asymme- 
tries, as expected. The systematic uncertainties in tag- 
ging (cr Q = 0.004, a c = 0.003) are estimated allowing 
for different tagging efficiencies between B° and B a and 
for different At resolutions for correctly and incorrectly 
tagged events. We also account for uncertainties on the 
background (a a = 0.001, a c — 0.003) by varying the ef- 
fective lifetimes, dilutions, toes shape parameters and 
signal fractions, and background CP asymmetry up to 
five times the expected CP asymmetry for signal. 

The results can be interpreted in terms of sin(2/3 + 7) 
(Eq. 0|| if the decay amplitude ratios A*\ expected to 
be Y^ub^cdl^udYc^ ~ 0-02, are known. Such small am- 
plitude ratios cannot be determined from B°^D < --*^ ± t: t 
events directly, because the current data sample is too 



small. We estimate r 
lation r-*' — tan 9c 



using the SU(3) symmetry re- 
5]. From the 



8(B°->Di, 



measurements of the Cabibbo angle tan 9c = 0.2250 ± 
0.0027 H3, the branching fractions B{B°^D-i: + ) 
= (0.30 ± 0.04)% [13, B(B°->D*-n+) = (0.276 ± 
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FIG. 2: Distributions of At for the B°~ >D ± 7r =F (a-d) and _B°— >D* ± 7r :F (e-h) candidates tagged with leptons, split by B tagging 
flavor and reconstructed final state. The lines are fit projections and hatched regions represent background. 



0.021)% 0, B(B°~>D+ir-) = (2.7±^±0.8) x 10" 5 0, 
B(B°^D*+7r-) = (1.9±iJ ± 0.5) x 10" 5 [3 , and 
from calculations of the decay constant ratios JdJId 
= 1.11 ±0.01 and f D *JfD* = 1.10 ±0.02 we obtain 

r = 0.019 ±0.004, r* = 0.017lg$$ . (6) 

To obtain sin(2/3 + 7), we minimize the \ 2 

X 2 (2/3 + 7, r W) = £ (^7^) ' + A(r«) , (7) 

where = a, a*, ci ep , C[* ep are functions of the physics pa- 
rameters (Eq.0J|, and the corresponding measured 
values. A(r(*') is a continuous function that is set equal 
to within 30% of the estimated (Eq. and is 
an offset quadratic outside this range, with the errors in 
Eq.® The additional 30% error attributed on is due 
to the unknown theoretical uncertainty on the validity of 
the SU(3) symmetry assumption and to neglecting W- 
cxchange contributions to A{B°^D^ + tt^). This error 
estimate is consistent with the spread in r^*> obtained 
using a variety of theoretical models 0] . The <ii are the 
quadratic sums of the statistical and systematic uncer- 
tainties in Eq. Correlations between the x% , at most 
28%, have negligible influence on the results of this anal- 
ysis. The simultaneous analysis of two B decay modes 
allows one to extract | sin(2/3 + 7)|. 

Figure [3] shows the minimum \ 2 for each value of 
I sin(2/3±7)|. The absolute minimum occurs for | sin (2/3+ 
7)| = 0.98, where Xmin/ d - oi - = 0.44/1. The values of 
r^*' that minimize the % 2 are consistent with the in- 
put values within their statistical errors. Because of the 
large uncertainties on the fit parameters and their limited 
physical range, the x 2 curve is non-parabolic. Thus to ob- 
tain a probabilistic interpretation to the results, we con- 
sider, for each of many values of sin(2/3+7), a large num- 
ber of simulated experiments with the same characteris- 
tics as the data. We compute the consistency of the data 
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FIG. 3: Dependence of % 2 on I sin(2/3 + 7)! (top) and of the 
frequentist confidence level of the agreement of the data with 
expectations as a function of the hypothesis on | sin(2/3 + 7) 
(bottom). The assumptions on r and r* are contained in the 
definition of \ 2 (Eq. - The dashed horizontal lines indicate 
the 68% and 83% confidence levels (defined in the text). 



with a given value of sin(2/3±7) by counting the fraction 
of simulated experiments for which x 2 (sin(2/3+7)) — Xmin 
is smaller than it is in the data. This fraction, the fre- 
quentist confidence level, is shown in the lower portion of 
Fig. 01 from which we read that | sin(2/3 + 7)! > 0.69 at 
68% C.L. We exclude the hypothesis of no CP violation 
(sin(2/3 + 7) = 0) at 83% confidence level. In order to 
study the impact of the assumed theoretical error on r^*' , 
we doubled it to 60% and we found that the lower limit 
on I sin(2/3 + 7)] at 68% C.L. drops from 0.69 to 0.60. 

In conclusion, we studied the time-dependent 
CP-violating asymmetries in fully reconstructed 
B^^D 1 **^^ decays, and measured the CP-violating 
parameters listed in Eq. [S] With some theoretical as- 
sumptions, we interpret the result in terms of sin(2/3±7) 
and we find that | sin(2/3 + 7)| > 0.69 at 68% C.L. and 
that sin(2/3 + 7) = is excluded at 83% C.L. 
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